Energy Spectrum of Superfluid Turbulence without Normal Fluid 
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The energy spectrum of the superfluid turbulence without the normal fluid is studied numerically 
under the vortex filament model. Time evolution of the Taylor-Green vortex is calculated under the 
full nonlocal Biot-Savart law. It is shown that for k < 2n/l, the energy spectrum is very similar 
to the Kolmogorov's -5/3 law which is the most important statistical property of the conventional 
turbulence, where k is the wave number of the Fourier component of the velocity field and I the 
average intervortex spacing. The vortex length distribution becomes to obey a scaling property 
refiecting the self-similarity of the tangle. 
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Particular attention has been focused recently on the 
similarity between superfluid turbulence and conven- 
tional turbulence [1-3]. Early work of the superfluid tur- 
bulence has been concerned with the counterflow where 
the normal fluid and superfluid flow oppositely [4], hav- 
ing no classical analog. However Stalp et al. studied re- 
cently the superfluid turbulence produced by the towed 
grid, thus finding the similarity between the superfluid 
turbulence and the conventional turbulence above 1.4 K 
[2] . They observed indirectly the Kolmogorov law which 
is one of the most important statistical properties of the 
conventional turbulence. This is understood by the idea 
that the superfluid and the normal fluid are likely to be 
coupled together by the mutual friction between them, 
and to behave like a conventional fluid [1,5]. Since the 
normal fluid is negligible at mK temperatures, an impor- 
tant question arises: even free from the normal fluid, is 
the superfluid turbulence still similar to the conventional 
turbulence or not? 

As the physical model to describe the vortex dynamics 
in superfluid He at very low temperatures, two types of 
models are well-known: the Gross-Pitaevskii (GP) equa- 
tion which describes the motion of a weakly interacting 
Bose condensate, and the vortex filament model governed 
by the incompressible Euler dynamics. The former re- 
duces to the Euler vortex filament model when variations 
of the wave function over scales of the order of the su- 
perfluid healing length are neglected. The GP equation 
includes such complicated compressible effects as the ra- 
diation of sound from the vortex lines [6,7], the vortex- 
sound interactions, etc. In order to consider the intrinsic 
property of superfluid turbulence in a simpler situation, 
we study the energy spectrum of the 3D velocity field in- 
duced by the vortex tangle in the absence of the normal 
fluid under the vortex filament model. 

The energy spectrum of the vortices in superfluid was 
numerically calculated by other authors. Nore et al. 
studied the energy spectrum of the decaying superfluid 
turbulence by using the GP equation, and finding the 



transient spectrum for small k has the Kolmogorov law 
[3] . However at late stage some complicated compressible 
effects become dominant. On the other hand, an advan- 
tage of the vortex filament model compared with the GP 
equation is the foUowings. First this model enables us to 
calculate the energy spectrum free from such complicated 
effects. Secondly some physical quantities, e.g. a total 
line length of vortices and a vortex length distribution, 
can be calculated easily, so that the relation between the 
statistical property in the wave number space and the 
self-similarity of the vortex tangle in the real space can 
be discussed. Under the vortex filament model, the en- 
ergy spectra were reported on the Kelvin wave cascade 
[8,10], which was limited to a few vortices, not a dense 
tangle. Thus the calculation of the energy spectrum of a 
dense tangle under the vortex filament model is expected. 

For superfluid ^He, the vortex filament model is very 
useful, because the vortex core radius uq ~ 10~*cm is 
microscopic and the circulation k = 9.97 x 10~^cm^/sec 
is fixed by quantum constraint. Since the Helmholtz's 
theorem for a prefect fiuid states that the vortex moves 
with the velocity produced by themselves, the dynamics 
is governed by the Biot-Savart law [9] . The velocity field 
due to the Biot-Savart law is divided into two parts: one 
is the localized induction field determined by a local cur- 
vature of vortex line, and the other is the nonlocal field 
obtained by carrying out the integral of the Biot-Savart 
law along the rest of the filament. When the dynamics 
of a dense vortex tangle is calculated numerically, the 
nonlocal velocity field is usually neglected (the localized 
induction approximation) [11,9]. However the dynamics 
in this work is calculated by the full nonlocal Biot-Savart 
law, because the long range effects can be important on 
this problem. A vortex filament is represented by a sin- 
gle string of points at a distance apart. When two 
vortices approach within A^, it is assumed that they are 
reconnected. The computational sample is taken to be 
a cube of size D — 1.0 cm. This calculation assumes 
the walls of the cube to be smooth and takes account of 
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image vortices so that the boundary condition may be 
satisfied. This calculation of the dynamics is made by 
the resolution = 1.83 x lO'^ cm and M = 4.0 x 10"^ 
sec. 

The energy spectrum E{k) is defined as E = 
Jg°° dkE{k), where E is the kinetic energy per unit mass 
and k is the wave number of the velocity field. In our 
previous papers we derived the energy spectrum under 
the vortex filament model [12]: 

xs'(a)-s'(6)e"''='^^^^^^"*^^=^\ (1) 

where dflk denotes the volume element k^ sin dkddkd(f>k 
in spherical coordinates. Here a vortex filament is repre- 
sented by the parametric form s = i), where s refers 
to a point on the filament, the prime denotes differentia- 
tion with respect to the arc length ^ and the integration 
is taken along the filament. The energy spectrum E(k) 
is calculated for the vortex configuration s(^) obtained 
by the simulation of the dynamics. 

Figure 1 shows the decay of the vortex tangle with- 
out the mutual friction [13]. As the initial configuration 
of vortices, we use the Taylor-Green vortex [3] (Fig. 1 
(a)). These initial vortices are highly polarized. How- 
ever through the chaotic dynamics which includes lots of 
reconnections, the vortices become a homogeneous and 
isotropic vortex tangle (Fig. 1 (b), (c) and (d)). 




(c) (d) 
FIG. 1. Time evolution of the vortex tangle at t=0 sec(a), 
t=3Q.O sec(b), t=50.0 sec(c) and t=7Q.O sec(d). 

First we discuss the transient behavior of the k depen- 
dence of the energy spectrum E{k). The energy spectra 
calculated from each configuration in Fig. 1 are shown 
in Fig. 2. It is shown that the slope is changed about at 
k = 2^/1, where I is the average intervortex spacing. The 



energy spectrum for k > 2Tr/l has k~^ behavior which 
comes from the velocity field near each vortex line [3,14], 
though the random vortices compose the turbulent ve- 
locity field. On the contrary, the spectrum for k < 2tt/1 
is strongly affected by the random vortex configuration. 
At t = sec, the spectrum has a large peak at the small- 
est wave number, being fiat in the intermediate range 
because there are only large vortices and no short-scale 
structure on them. Figure 2 shows that as the vortices 
become the homogeneous and isotropic vortex tangle the 
slope for k < 2tt/1 converges to the Kolmogorov form 
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FIG. 2. The energy spectra of the vortex tangle at t=0 sec 
(dashed), i=30.0 sec (dot-dashed), t=50.0 sec (long-dashed) 
and t=70.0 sec (solid). 

The Kolmogorov law can be derived from the argument 
based on the picture of the cascade process [15] . In the in- 
ertial range the kinetic energy is transferred steadily from 
small k to large k without dissipation, and dissipated at 
the end of the inertial range. Thus, for the steady state, 
the energy dissipation rate e = ^dE/dt can be identified 
with the energy flux in the inertial range. Then the en- 
ergy spectrum depends only on the wave number k and 
the energy dissipation rate e = —dE/dt, which leads to 
the Kolmogorov spectrum E{k) — Ce^/'^fc^^/'^. Here C is 
the (dimensionless) Kolmogorov constant of order unity. 

In our previous papers, the cascade process without 
the mutual friction in superfluid turbulence was dis- 
cussed [9]. Through lots of reconnections, the vortex 
tangle breaks up to smaller ones and this process pro- 
ceeds self- similarly, and in our calculation the smallest 
vortex whose size is the order of the numerical space res- 
olution is eliminated by the cutoff procedure. This 
resolution in our calculation is much larger than the dissi- 
pative scale in real system. However this numerical cutoff 
can be justified, because the cascade process at a small 
scale proceeds much faster than that at a large scale. Ac- 
tually we showed the decay rate of the density of vortices 
was almost independent of the cutoff scale A^. Figure 3 
shows the energy dissipation rate e due to the cutoff pro- 
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cedure in the dynamics of Fig. 1. After 70 sec, the tangle 
becomes isotropic and homogeneous losing the memory 
of the initial configuration, so the change of the dissipa- 
tion rate becomes slow free from the artifact of the early 
large dissipation. 
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FIG. 3. The energy dissipation rate e which calculated from 
the dynamics of Fig. 1. 

Next we compare quantitatively the energy spectrum 
at 70 sec with the Kolmogorov law E{k) = Ce^/^k~^/^. 
The Kolmogorov constant C is known as the parameter 
of order unity. Here we use C = 1 and e = 1.287 x 10~^ 
cm^/sec^. Then we can determine uniquely the energy 
spectrum. Figure 4 shows the comparison between the 
energy spectrum at f = 70 sec and the Kolmogorov law 
with C = 1 and e = 1.287 x 10"^ cm^/sec^. The en- 
ergy spectrum for k < 2tt/1 is consistent with the Kol- 
mogorov law not only on the wave number dependence 
but also on the absolute value. The dissipative mech- 
anism due to the cutoff works only at the largest wave 
number k ^ 27r/A^ = 343 cm^^ . However the energy 
spectrum at small k region is determined by that dissi- 
pation rate. This result supports just the picture of the 
inertial range. Although our spectrum has k^^ region 
between the Kolmogorov region and the dissipative wave 
number, the energy flux exists also in this region. The 
spectrum in this region includes the contribution coming 
from each vortex line and that of the energy flux, while 
the former is dominant [14]. The cascade process in this 
region will be discussed elsewhere. 

The Kolmogorov law is the scaling property in k space, 
and closely related with the self-similarity of the turbu- 
lent velocity field in the real space. We devote the rest 
part of this paper to the following question: is this scal- 
ing property in k space related with the self-similarity 
of the vortex tangle in real space or not? However, in 
conventional turbulence it is very difficult to discuss this 
problem, because the viscous diffusion of vorticity makes 
the vortex configuration obscure. On the contrary, the 
characters of the superfluid turbulence are the definite- 



ness of the vortex line due to the absence of the viscosity 
and the fixed circulation by the quantum effect. These 
characters allow us to describe the system by the topo- 
logical configuration of a vortex tangle. Hence, in order 
to discuss the self-similarity in a real space, it is mean- 
ingful to investigate a vortex length distribution (VLD) 
where n{x)Ax represents the number of the vor- 
tices whose length is from a; to a; -|- Ax. 
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FIG. 4. Comparison of the energy spectrum (solid line) 
at t = 70 sec with the Kolmogorov law E{k) = Ce^^'^k''^^^ 
(dotted line) with C = 1 and e = 1.287 x lO"*^ cm^/sec^. 

In order to suppress the fluctuation of the VLD, we 
use the time averaging procedure in the interval 4.0 sec. 
Figure 5 shows the averaged VLD at each time, where 
the largest scale in x axis is the size of the cube and the 
smallest scale is the length of the numerical space res- 
olution. At large scale the VLD is strongly affected by 
the effect of the boundary, and at small scale decreases 
rapidly by the effect of the cutoff procedure. As the vor- 
tices approach the homogeneous and isotropic tangle, in 
the intermediate range the VLD becomes to obey a scal- 
ing property n{x) oc x". Using the least squares flts, we 
determine this scaling exponent a at t = 60.0 sec in the 
region of 0.15 < a; < 0.65 [cm], and find a = — 1.34±0.18. 

In conclusion we have studied numerically the energy 
spectrum of the superfluid turbulence without the mutual 
friction. For k > 27r//, the spectrum can be attributed to 
the contribution of each vortex line. For k < 2tt/1, as the 
vortices approach the homogeneous and isotropic tangle, 
the slope of the spectrum converges to the Kolmogorov 
form /c"^/^. By using the energy dissipation rate due to 
the elimination of the smallest vortices, the spectrum for 

< 27r/Z is consistent with the Kolmogorov law not only 
on the wave number dependence but also on the absolute 
value. The VLD n{x) becomes to obey a scaling property 
reflecting the self-similarity of the tangle. 
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FIG. 5. The VLD n{x)Ax at t=0 sec (dashed), i=20.0 sec 
(loug-dashed), t=40.0 sec (dot-dashed) and t=60.0 sec (solid) 
for Aa;=0.05. The dashed line is determined by the least 
squares fits aX t = 60.0 sec. 
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[14] The energy spectrum of a straight vortex line can bo de- 
rived. Since the velocity field of a straight vortex lino is 
it/2nr, the kinetic energy of a vortex line per unit mass is 
calculated as E = (1/2) / v^{r)d^r = {K^L/4n)ln{D/a), 
where L is the total line length per unit volume and r is 
the distance from the vortex lino. Since the energy spec- 
trum E{k) = Ck" can be defined as E = J^^j^ E{k)dk, 
the parameters a and C are determined as a = — 1 and 
C — k^L/Atv. For example, at t =70 sec, L = 151.9 cm~'^, 
so that E{k)iinc = (K^L/47r)fc"^ = 1.20 x 10"^ cm'Vsec^ 
for k —100 cm^^. This is consistent with our numerical 
result E{k) = 1.07 x lO"'' cm^/sec^ for k =100 cm"^ in 
Fig. 4. 

[15] Sec, for example, U. Frisch, Turbulence (Cambridge Uni- 
versity Press, Cambridge, England, 1995). 
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